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Abstract 

We study self-similar solutions of NLS-type dynamical systems. Lagrangian approach 
is used to show that they can be reduced to three canonical forms, which are related by 
■ Miura transformations. The fourth Painleve equation (PIV) is central in our consideration 

' - it connects Heisenberg model, Volterra model and Toda model to each other. The 

^\ , connection between the rational solutions of PIV and Coulomb gas in a parabolic potential 

' is established. We discuss also the possibility to obtain an exact solution for optical soliton 

On ' i.e. of the NLS equation with time-dependent dispersion. 

ON ■ 

1 Introduction 

The self-similar solutions are a very useful instrument to investigate some important physical 
Q , applications of integrable systems. 

^ ' There are many examples of these applications. 

The importance of this consideration follows from the fact that the asymptotic form of 
any solution is a self-similar one. Because we deal with integrable systems we can build exact 
self-similar solutions from their asymptotic. This approach looks like a form of Painleve test. 
' In this article we establish a connection between the self-similar solutions of the NLS-type 

dynamical systems and the fourth Painleve equation (PIV). 

The purpose of this article is to obtain the equations for the self-similar solutions of NLS- 
type equations. We show that any member of this NLS family can be reduced by translation 
and dilation transformations to three fundamental canonical forms. Miura transformations 
among these canonical forms and Backlund transformations will be used to construct new 
solutions. We show that the Miura transformations are preserved also under the self-similar 
reduction. 

We find also an example of an exact self-similar solution for the NLS equation with time- 
dependent dispersion. 

The structure of the article is the following. In the next section we consider the Hamil- 
tonian formalism for a NLS-like family of equations and we present their possible canonical 
forms and Miura transformations relating them. In Section 3 we derive the equations for the 
corresponding self-similar solutions and show that all possible cases are connected with the 
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Painleve IV equation (PIV). We present a few exact solutions, interesting for some physi- 
cal applications, namely a Toda system corresponding to a "collapsed" solution of the NLS 
equations and the Heisenberg system describing spin dynamics. 

2 Canonical forms and Miura transformations 

Hamiltonians of NLS-type dynamical systems have the form: 

H = PxQx+ eP^Qx + c^PQl + PQl + iP^Qx + 5P\ (1) 

where e, a, f3, 7 and S are arbitrary parameters. 
By using translations 

P-^p-po, Q ^ q - xqox -tqot, X~^x~vt, (2) 

and dilations of the canonical variables P, Q and space X and by re-scaling the Hamiltonian, 
the Hamiltonian (|l|) can be transformed into one of the following canonical forms. 

2.1 Canonical forms 

I. Heisenberg model e 7^ Hamiltonian (|l|) by means of the transformation 

2e 2e e e 

reduces to 

2 2 

H = p^q., -f ep\l + {5- + (/5 - (3) 

p 1 

By the dilations p — *■ and iJ — > — H and by renaming the constants we have 

e e 

H ^p^q^-p^ql + 6p^ + I3ql. (4) 

The generic case with 5^0 and /3 7^ is considered in Section 3.3. Of special interest 
are the two special following cases. 

A Isotropic Heisenberg model (/? = 0, (5 = 0), i.e. 

H = p^q^-p'^ql. (5) 

B Anisotropic Heisenberg model [(i ^ Q, 5 — 0), i.e. 

H ^p,q,-p^ql + Pql. (6) 

The other case {S ^ 0, j3 = 0) is obtained from this one by the transformation 
p^ qx, (3^5. 
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II. Volterra model (e = 0, a ^ 0, 77^ 0.) Hamiltonian reduces to 

H ^p^q^+pq^+p^q^. (7) 

by the transformation 

p = p-tL^ Q = q^-x t--—t, X^x + 2[— + ^]t (8) 

a 7 7^ a \ 1 J 



and by the dilation 



p q 
p^ -, g -> -. 

7 a 



III. Toda model (e = 0, a ^ 0, 7 = 0.) Hamiltonian (|^) reduces to 

H ^p,q,+pql+p^. (9) 
by the sequence of transformations 

P = p--, Q = q~2^t (10) 
a a 

P^T, 9^-- (11) 
a 

Note that p ^ qx <S=^ a ^ 7, /3 ^ (5 . Other cases lead to linear problems. 
2.2 Miura transformations 

One can easily verify that the NLS-type systems are connected by Miura transformations 
according to the following scheme: 

isotropic Heisenberg model (I) Volterra model (II) Toda model (III) (12) 

Px 

I ^ II : -pqx-\ > P, -pQx ^ Qx 

P 

II ^ III : pqx +Px P, P + Qx-^Qx- 

3 Self-similar solutions 

Let us choose solutions for the NLS-type systems in the form 

qx^b{t)y{0, p^c{t)z{0, ^ = xa{t), (13) 

To obtain equations for y and z one can consider p as a lagrangian coordinate, ^ as the 
new time, substitute ansatz ([l3| ) into the Lagrangian 

Li=pqt-H, (14) 

and require that the dependence on t factorizes (i.e. the dependence on t of all terms in 
Lagrangian must be the same). Then, by considering the variation of the corresponding 
action with respect to z one gets an equation for y and z. 
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In the same way one considers g as a lagrangian coordinate and requires that the t depen- 
dence in the Lagrangian 

L2 - -qpt - H, (15) 

factorizes. Finally, evaluating the variation of the action with respect to y one obtains the 
second equation for y and z. 

Lagrangians (^4|) and ( p^ ) can be rewritten as 

Li — Loi — V, 

where the "bare" parts of the Lagrangians 

Loi = pqt - PxQx, Lq2 = -qpt - PxQx (16) 
are the same in all models. Requiring the factorization of time for both Loi one gets 

a' = 2a^ 5 = a^ c = a^, (17) 
and then the "bare" parts of the actions have the form 

^01 ^ Jdt b{t) c{t) J rfe[2(^ - l)zY + 2^zy ~ z'y], Y'^y (18) 

S02 ^ J ^(0 "=(0 J d^[-2XzY - 2^z'Y - z'y]. (19) 

By variation on z of the first action and on Y of the second action one derives the system of 
equations 

y' + 2^y + 2{^-l)Y + c,^^ 

-z' + 2e^ + 2(A-l)Z + C2 = %^, Z' = z ^^"^ 

for and A such that V{p,qx) — V{cz,by) — abcv{z,y). Note that, respectively, for v ^ 1 
and for A 7^ 1 the constant ci and C2 can be omitted. Note also that the system is lagrangian 
only for A + = 2, since only in this case the two actions in dl^ ) and (^9|) differ by a total 
derivative with respect to ^. The next considerations depend on the model (i.e. on V). 

4 Volterra model 

In this case 

V=pql+ p\x ^i^ = X = l (21) 
and we obtain the Lagrangian system {L = y' z + 2S^yz — y^z — z^y + c\z + C2j/) 

y' + 2iy + ci = 2/2 _,_ ^zy 
-z' + 2iz + C2 = 2yz + 

which is an equivalent form of the PIV ODE . In fact inserting z in L from the first equation 
of ( p2| ) and writing 

w(x,a,6) = ?/(-^,ci,C2), a=ici-C2 + l, b = -^cj (23) 
we obtain the PIV Lagrangian 

L= — +w^+ AuPx + Aw{x^ - a) - — (24) 
w w 
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4.1 Coulomb gas. Rational solutions of PIV. 

Another useful representation for PIV is 

5i = 51(53 - 52) + ai, 52 = 52(51 - 53) + "2, 53 = 53(52 " 51) + "3, 
with additional first integral 

51 + 52 + 53 = x^, 7 = ai + a2 + 03. 
It easy to see that all of 5^ = j = 1, 2, 3 are solutions of PIV if 

c. ct 
7 = -2, " = ^= — ^' Cj=aj+2-ai+i 



It is well known (see ||3[ for example) that system (g5|) has the rational solutions if 
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-7j, 7j e Z, 71 + 72 + 73 = 3, 71 = 72 = 73 ( mod 3). 



^ 3 
We write 

5i = -^73^ + log(Qj+i/Qj-i) 

where - polynomials and 

(g;' - xq; + m,Q,)g,+i + (g^+i + xQ'^+, - m,+iQ,+^)Q, = 2g;g;+i, (7 = 3). 

For example 

a-(-l,2,2), g = (x,l,l), a = (-l,5,-l), g = (x^ - 1, 1, + 1), 

a = (1,4,-2), g= (a;,l,a;2 + l), a =(5, -4, 2), g = (a;, x'' + 2a;2 - 1, a;^ + 1). 
It is important for us that any divergence of PIV has a form 



c 



- Co, = 1 



It is well known also that the asymptotical behavior of functions gi defines by 

51 -2C, 52, 53 ^ ^ ^ oo, 
We are looking for rational solutions of the PIV equation in a form 

where 5 = | in a case ( |3^ ) and 6 — oi S — 2 in (|33|). 



Substituting (!§ to (|l]) one has 

iS~m^J2T^^ c, = ±l. (35) 

Equation (^5|) describes the gas of charged particles with Coulomb interaction in a parabolic 
potential, and the level statistics in a hermitian case because one can obtain ( |35| ) by variation 
of a functional 

^ = E ^(1 - ^)^'^' + E ^^^^ - ^^■)- (36) 

In the case when all charges equal Ci = 1, (or Ci = — 1, ) one have a famous solution of PIV 

»^-2? + f|, (37) 

where Hn{^)— Hermitian polynomials (c.,[^). So, the roots of i7„ coincide with distribution 
of the positive (negative) charges in the potential (|3^). 

For example, solution (5(1,4,2) (see (|30|)) defines three solutions of ( ^5|) if 5 = |, two 
are "trivial" (Hermitian polynomials) and third describes a system of two positive and one 
negative charges: 

ei,2 = ±i, 6 = 0, Ci,2 = l, C3 = -l. (38) 



In general case one have three solutions of (^5|): all positive charges coincide with the roots 
of Qj+i, all negative charges coincide with the roots of Qj-i; numbers of charges are N± = 
deg Qj±i. 



5 Heisenberg model 

In this case 

V = -p\l ^X^l-y, (39) 



and from (20) we have 

y' + 2^y + 2{v^l)Y = -2zy^ 

-z' + 2^z - 2uZ ^ -2z^y. ^ ' 

This system admits the conservation law 

zy' - yz' + 2^zy + 3zY = (41) 

where /i is an arbitrary constant. Inserting z from the first equation of (^) into this conser- 
vation law we get a third order ODE satisfied by Y 



-— + 8(:/- 1)C— + 6 i^- 1) — 

y 2 y^ y y^ 



+ 8{iy-l)^-+6{iy-iy^ + 2C^ + 2{v-ii)^ 0. (42) 
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By setting 

Y = el'd^'wii')^ (43) 

we have 



^ ..2 I o^^.. T \ ^ I e/.. 1 \2 1 , 



-— - -w' + 8e(i. - 1)- + 6(j/ - ly— + 2e + 2i,y-^i)^ 0. (44) 

^ ■■ I ^ / T/i Td-^ 



i« 2 2 w 

Finally, by means of the transformation 

one gets the PIV ODE 



(45) 



/" = + If + 4?/^ + 2{e + ^-n)f- 2^^y^- (46) 

5.1 Self-similar solutions for spin dynamic 

If we introduce the vector ~^ 

Si^p{q^-l) + q (47) 

52 = tp{q^ + I) + tq (48) 

53 ^2pq+l (49) 

of unit length, i.e. — 1, the Hamiltonian (||) for the isotropic Heisenberg model can be 
rewritten as 

H = -\'^l (50) 

and the corresponding evolution equations 

Pt = {Px - ^p'^Qx)^ (51) 

qt = -qxx- 2qlp. (52) 

as 

f§ t = ^ X XX (53) 

which up to a change to the physical time t it is the system of equations describing the 
Heisenberg Ferromagnetic model. 

In the previous section we show that the general self-similar solution is related to the PIV 
Painlevc transcendents. It is interesting to note that the requirement 

5-3 = cr, a = ±1 (54) 

is compatible with the existence of self-similar solutions. More precisely, the solutions are 
related to the parabolic cylindric functions. In fact from 

53 = 1 + 2zY (55) 
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inserting z from the first equation of the system (|4C|) and using ( [4 3D and (|45|) we have 

/' = 2a{v - 1) + 2U + /' (56) 
and one can verify directly that these / are solutions of the PIV equation (^) for 

v-y^^\j^a{v-\). (57) 
The Riccati equation ( ^6|) can be linearized by writing 

/ = (58) 

The -0 satisfy the Weber equation 

V-" - - 2a{v - 1) - l) i/) = 0, (59) 

whose general solution is a linear combination of the parabolic cylindric functions 

^1 -^.(.-i)(\/2e) (60) 
^2 = ^.(.-i)(-V2C). (61) 

5.2 Deformation of Heisenberg model 

First, let us note that the evolution equations corresponding to the general Hamiltonian (Q) 
with /3 7^ and S ^0 do not admit self-similar solutions of the form (]l3|). In the case 

V = ~p\l + Pql (62) 

the evolution equations are 

Pt = [Px - Sp^gx + '2Pqx)x (63) 

qt = -qxx - 2qlp (64) 

and self-similar solutions are obtained for v — 1, A = 0. 
We have 

2z = g' - 2^g + ag^ (65) 

y = I (66) 

and 

g'" = 4^5 + icjg'g' - ^c.^gg' - ic,g' + (^^^^) (67) 

with ci a constant of integration. By introducing g — we have 

w'" = ~w'w" + 4^ + 4/3w'e-2'" + 

+ iclw'e^"" - Aci£,w'e^ - 4cie™ (68) 
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that can be easily integrated once to 



w" = + 2e - 2/3e-2- + 2 + 

+ ^c^e^"' -4ciCe"'. (69) 
where 7 is a constant of integration. Coming back to g we get 

25 9 
+ ^c?5^-W. (70) 

For ci 7^ changing ci.g g we get the PIV ODE equation. Dividing ( |70| ) by differentiating 
with respect to ^ and then by inserting in the obtained equation g" from ( |70| ) one gets 

5"' = + 4 (e^ + 7) 5' + ecjg^g' _ Uci^gg' - ^c^g^ (71) 

which for ci — 1, together with (|6^), can be considered an equivalent form of the PIV equation. 
Of special interest is the case Ci = 0, i.e. 

g'" ^4^g + i{e+7)9'- (72) 
One can easily verify that its general solution is given by 

9 = diipf + d2^2 + dz'ipiil)2 (73) 

where dj {j = 1, 2, 3) are arbitrary constants and ipi, tp2 are two independent solutions of the 
Weber equation 

r-{e+i)i^^o. (74) 

e.g. the parabolic cylindric functions 

= D_2+i(\/20 (75) 

^1^2^ D_J + l{-^/2C}. (76) 

5.3 Landau-Lifshitz model. 

The Hamiltonian of the Landau-Lifshitz model 

H=p^q^-p\gl+r{q))-p'-^~'-^, = (77) 

corresponds to the "diffusion" spin dynamics 



iSt^[Sx{S.,, + JS)], (78) 



where we introduced the "canonical" parametrization of the vector S 

S = ^2, ^2) = {p{q^ -l) + q, ^P{q^ + 1) + ^q, 2pq + 1), r{q) - -^iSpSSp). (79) 
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One can obtain another parametrization of S from a spectral problem of the Toda model 



(see (|87|)): 

w{x)S3 ^ {ipiil:2)x, w{x){Si+iS2) = i{ipi)x, w{x)[Si - 182) ^ i[ipl)x- (80) 

where w = - i'2,x'ipi- 

It is important that the Frenet equations 

Sx^kN, N,^-kS + xB, R,^-xN, {Sj)^l, (81) 

where B — S x N— binormal, can be obtained from ( p7| ) if 

fc2 = 4p, ix + qx + (log k)^ = A, (82) 

6 Toda model 

In this case 



V ^pql+p^ ^ X^2, ^ = 1, (83) 



From (EOl) we have 



y' + 2^y + C = + 2z 
-z' + + 2Z = 2zy. 

Excluding y one obtains the third-order equation on Z 



(84) 



Z'" = — ^ + 2{Z'i^ j)- 4Z'2 + 2(C + 2)Z' (85) 

2Z Z 

which is integrable to 

Z"2+4Z'3-4(C + 2)Z'2+4i:)Z' = 4(eZ'-Z)2, Z? = C2 (ci + 2) . (86) 
Zero-curvature representation for Toda model (^) has a form 

i^xx + (gx - A)V'a; -I- p^A = 0, V't = (fe - A)?/;^; - pV' (87) 
One can substitute self-similar form ( ]8^ ) to obtain zero-curvature representation for PIV: 

^« + (j/ - A)Vc + zv- = 0, AV-A = (A + e + y)V'« - zV-, (88) 

where w = y + ^ is solution of PIV. 

As we shall show in the following section, (^) is another equivalent form of the PIV 
equation. It is interesting to note that it can be considered a conservation law of the system 
(p4[). In fact inserting Z from the second equation of (|8^) into (|6|) we get 

z\y - 20 + z{y - 2^^ ~ + {C + 2)z = D. (89) 

The function Z can be also expressed in terms of three PIV transcendent yniS,) [n G -^3) 
according to the formula 

2Z = yiy2y3 - "12/2 + ^22/1 (90) 



10 



where the parameters are defined by 

C — a2 — ai — 2, D = —aia2- (91) 
and the functions Hn are related by the equation 

yi + y2 + 2/3 = 2^ (92) 

and satisfy the PIV equations 

y'n-^jlvl- ^^Vl + + ln)yn + (93) 



with 



7n = ""^^ : ""-^ (94) 



Pn = -al (95) 

ai + a2 + = —2. (96) 

This result can be obtained using the connection among nonlinear chains and Painleve equa- 
tions found in j^. In fact let us consider the nonlinear chain 

y'n + VniVn+l " Vn-l) + Q!„ = 0, Tl £ Z3 (97) 

where the constants an are normalized as follows 

ai + (12 + as = -2 (98) 



and the y„'s satisfy (92) 



Easy computation shows that indeed the y^s satisfy the PIV equations (|9^). Let us now 
introduce 



Then we have 



and 



h = yiy2 (99) 
Z' = h (100) 



2/ 



h 

From these equations, eliminating yi and y2, we get finally that / satisfies 
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The system (^3|) is equivalent (up to the change t — > it) to a system of two coupled NLS 
equations. In fact if we rewrite it as 

Pt = {Px + '2pqx)x (103) 

qt = -Qxx + ql + 2p (104) 

and introduce 

u = e-« (105) 

p = uv (106) 

we get 

ut + u^x + 2u^v = (107) 

vt - v^^ - 2v'^u = 0. (108) 

7 NLS equation with time-dependent dispersion. 

7.1 Integrable case 

The special case of the NLS equation deserves a more detailed analysis and an extension of the 
study to some nonintegrable variants. Let us consider the NLS equation with time dependent 
coefficients 



tqt + f{t)qxx+g{t)\qW = Q. (109) 



Renaming the time 



and the function g(t) 



it can be rewritten as 



t 



f{t')dt' (110) 



m - 4i (111) 



fit) 



iqt + qxx+9{t)\q\' q^Q. (112) 

It admits self-similar solutions for g{t) of the form 

g{t) = a(4t)^ (113) 

where m € R and the constant coefficient a can be normalized to ±1, the sign being most 
important. The self-similar solution (up to a trivial rescaling and t > 0) is given by 

g=|g|e^«("'*), \q\' ^ 4k{tr+^Y'{fj), r^^k{t)x, k{t) ^ (114) 
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where 



(115) 



t 

2 

and Y satisfies the ODE 



0{v)^'^+mj ^dn, (116) 



2Y"'Y' - + Arj^Y'"^ - Arri^Y^ + 16aY'^ - 166*0 ^'^ = 0. (117) 

A Painleve analysis of this equations shows that it has no movable critical poles only for = 1 
and we recover equation (^) up to a trivial rescaling in order to pass from the variable ^ to 
the physical variable rj. In this case the integration factor Y" /Y''^ allows to integrate and to 
reduce it to the second order ODE 

Y^^ + 4 {rjY, - Y f + 2aY^{2Y, ~ aOn f + 2a^i% = 0. (118) 

with /i^ a constant of integration. This self-similar solution of the NLS equation was first 
obtained in Q in the case 6o = 0. Its properties can be studied following the same lines as in 

i- 

The solution Y can be expressed as 

-aY =^WiW~ + ^ [W^ - 2W, - /i^ + 1] + {W ~ t^) . (119) 



where W satisfies 



which up to a trivial chan ge o f variable is the PIV ODE. Equation (|119| ) can be considered 
an algebraic resolution of ( 118 ) which is quadratic in YJ,^. 
We have also 

-2«y„ = (i^-,)V^^^^4^. (^2^) 

Therefore for a — ~1 and for /i and W real the condition Y,, > is automatically satisfied. 
On the contrary for a = 1 and /x and W real we have 1^, < and we have not self-similar 
solutions. Let us then choose from now on a = — 1. 

We can also express W in terms of Y as follows 

^ ^ ^f^Yyv + {2Yr, + eo) {2Y + 7?^o) + Tjfi^ 
{2Y„ + eof+fi^ 

Starting from a suggestion by Bureau ^, pag. 210] one can solve algebraically ( |118| ) in an 
alternative way and, precisely, one can express Y in terms of another function W as follows 

^ ^ b 

4Y = hr! + -^ - 2W^ - AnW^ + 2 (ri^ ~ a)W + — (123) 
2W ^ ' AW 
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with W satisfying the equation 

WWr,r, = - QW'^ + SrjW^ -2{if-a- cr^) W"^ 

where 

and a, h and h are defined by 

+ 3(3/x^ - 6*^)0 + 6»o (Q/x^ + 6*0) = 

h=\{a- 200) . 
The roots of the algebraic equation for a are 

ai = -9q 

1 3 
012,3 = ± 2*M- 

and then 

61 = -2^2 



1 2 



/ii = -200 
^2,3 = -6*0 ± ijJ- 

W can be expressed in terms of Y. 
In the case ai = —9o one gets 

Wi = -rj--—-—^ 

and in the case 02,3 = ± |«M 

- ^ y^^ 



W2- 



2Yr, + eo±ilJL 2 2Yr,+0o±in 
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The W are related by the equation 

C2W2 



C3W3 



W2 ~Wi W3- Wi 



C3W3 - C2W2 



where 



C2,3 = - (6'o ± i/i) 



(136) 



(137) 



The Backlund transformation for the NLS equation induces a Backlund transformation for its 
self-similar solution Y. Precisely one gets a new self-similar solution Y which satisfies ( |118| ) 
with —t 2 — pu and which is given by 



y = y + i 



(Ai-1) (2r^ + 0o) +m' 



2 [-/iK,,, + {TYr, + 0o) (2y + r;0o) + m^"" 



(138) 



7.2 NLS with time-dependent dispersion. A non integrable case 

Let us consider the Lagrangian for NLS with a time-dependent dispersion 



L ~ dt dx 



i^>*-d(t)^:Vx + |(^*V')' 

The canonical transformation {tjj^^p*) {p,<f) where -0 = p^/^ e*"^ gives us 



dt / dx 



L 



The equations of motion are given by 



+ 2ddx{p<px) = 0, 

4>t +d(j)^ - -I- - — 

2 p 4 



.gp = 0. 



The continuity equation (141) is identically satisfied for p of the form 
P = A(t)/(0, cj, = - 



with A, /, and (j)Q arbitrary functions. 
The Lagrangian (|140|) transforms to 



AdX{ty 



(139) 

(140) 

(141) 
(142) 

(143) 



'U^\-dxf^ 



A KMX J 



V2dA2 



(144) 



The equations of motions are compatible with a dependence of / only on ^ for a time dependent 
dispersion d and function 0o satisfying the requirements 



Ad = a, ^ 
A 



-/3, d"d : 



(145) 
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with a and /3 constants. Then the equation of motion for / is given by 



f 2/2 



or in terms oi = f 



y" + -ev + -y^ + /3y = 0. 

a a 



The equation for the dispersion d can be integrated once to 

[d'Y = ed^ + W^. 



(146) 



(147) 



(148) 



If a'y > 0, e < we have periodical solutions for d{t), but the equation for / has no finite 
solutions; if 07 < 0, e > 0, one can obtain the solution 



4a^7 



cosh {t\/e). 



(149) 



We say we have obtain an exact self-similar solution of NLS - equation with time dependent 



dispersion because solutions of (147) have no divergences and have well-defined asymptotical 



behavior (Schrodinger equation for oscillator with cubic nonlinearity) . It is important that we 
have no exact solutions in a case of periodic dispersion. So, an optical soliton [p^ , [ll|] in the 
periodic case may be non-stable. 

We have obtain the self-similar equations for NLS-type dynamical systems. It is possible to 
apply the Lagrangian method to the integrable two-dimensional systems or to non-integrable 
systems. 



Appendices 

A Miura transformations for self-similar solutions 

The sequence of Miura transformations (^2|) leads to analogous ones for the self-similar solu- 
tions: 



isotropic Heisenberg model (I) Volterra model (II) Toda model (III) 



(150) 



I ^ II : 



z 

~zy H 



z, -zy^y, 2(1 - I/) = C2 - ci, n = Ci 



II III 



zy + z 



z, z + y^y, C = ci+C2, D ^ C2{ci + 2) 



For obtaining the relation among constants of different reduced systems, in the case / II, 
one inserts the Miura transformed solutions z and y into the system // getting for the first 
equation the conservation law (|4^) containing the constant ci and for the second equation the 
same conservation law containing the constant C2, while, in the case // —^ III, one inserts the 
Miura transformed solutions z and y into the conservation law (p9|) of system /// and then 
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uses the second equation of the system II and the PIV equation for the self-similar solution 
y of the system II. 

It is interesting to note that the Miura transformation II — > III by using the first equation 
of the system II can be inverted getting 

III ^ II ■■ y > z, — > y. 

B Lagrangians for Painleve equations 

The six Painleve equations (a, 6, c, d constants) 

PI: W:,^ = &w^ + x, (151) 
P2 : = 2w'^ + xw + a, (152) 

P3 : w^^ = ^- — + + cw^ + -, 153 

w X X w 

PA: w^^ = -^ + — +Axw'^ + 2{x'^-a)w+-, (154) 
2w 2 w 



,1 1 . ^ Wx 2 , b. cw dwiw + l) 

— + -)wl - — + —{w - if {aw + - + — + ^ 

"2w w-l' X x^ ' ^ ^ w' x w-l 



PQ - w^^ = + ^— + ^—)wl - {- + ^— + ^—)w^ + (156) 
2w w — 1 IV — x x X — 1 w — X 



^ w{w — l){w — a;) 

x'^(x — 1)^ 



bx c , , , dx(x — 1) 

« + -3 + 7z— rT2(^-i) + 



w 



{w — 1)^ {w — xY 



can be derived by variation of the following Lagrangians: 

-^+2w^ + xw, (157) 
wl+w'^ + xw^ + 2aw, (158) 

-^x + aw-- + -{cw^-—), 159 
2w^ w 2 w' 

y,2 2b 

-^+w^+ Aw^x + Aw(x'^ -a) , (160) 

w w 

wl 4 b 2c 2dw 

^x+-{aw ) —tX, (161) 



LI : 


L 


L2 : 


L 


L3 : 


L 


LA : 


L 


L5 : 


L 


L6: 


L 



w{w — 1)^ ' x'^ w' w — 1 (u> — 1)^ 

x{x — 1) aw b c d 



2 w{w — l){w — x) x{x — l) w{x — l) x{w — l) w — x 
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